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Abstract 

We introduce an extended version of the fractional Ornstein-Uhlenbeck (FOU) process 
where the integrand is replaced by the exponential of an independent Levy process. We call 
the process the generalized fractional Ornstein-Uhlenbeck (GFOU) process. Alternatively, 
the process can be constructed from a generalized Ornstein-Uhlenbeck (GOU) process using 
an independent fractional Brownian motion (FBM) as integrator. We show that the GFOU 
process is well-defined by checking the existence of the integral included in the process, and 
investigate its properties. It is proved that the process has a stationary version and exhibits 
long memory. We also find that the process satisfies a certain stochastic differential equation. 
Our underlying intention is to introduce long memory into the GOU process which has short 
memory without losing the possibility of jumps. Note that both FOU and GOU processes 
have found application in a variety of fields as useful alternatives to the Ornstein-Uhlenbeck 
(OU) process. 

Keywords: Ornstein-Uhlenbeck processes, Levy process, Stochastic integral, Long memory, 
Fractional Brownian motion. 
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The fractional Brownian motion (FBM) is one of the most popular processes for constructing 



$_j ■ long-range dependent stochastic processes with continuous path and its fields of applications are 



very wide. To name just a few, we see FBM models in the fields of tele communications , signa l 



processes, environmental models and economics. A recent refe rence i s e.g., Doukhan et al.l (120031 ). 
Statistical methods for FBM have also been studied (see e.g. iBeranl (J1994J)). 
We review the definition and name properties of FBM. 

Definition 1.1 Let < H < 1. A fractional Brownian motion B H := {Bf} te ^ is a centered 
Gaussian process with B^ = and Cov (Bf , Bf) = \ {\t\ 2H + \s\ 2H - \t - s\ 2H ) , t,s G M. 



* Correspondence to Department of Mathematics, Keio University, 3-14-1 Hiyoshi Kohoku-ku, Yokohama 223- 
8522, Japan. E-mail:mmuneya@gmail.com. This research is supported by JSPS Research Fellowships for Young 
Scientists. 



From the definition FBM has stationary increments and is self-similar with index H, i.e., for c > 
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{Bf}t 



= {c H B^} te ^ where = denotes equality of all finite dimensional distributions. While 
with H = 1/2 is a two-sided Brownian motion (BM) and has independent increments, 
with H G (0, =) U (|, 1] has dependent increments. For < h < s and tGK and N EN, 
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Thus, £> H with //" e (-, 1] has a long memory property, namely, ^^1 Thin h) = oo. Finally B H is 



known to have bounded p- variation for 1/i? < p < oo (see Prop osition 2.2 oflMikosch and Norvaisa 
( 2000)). For a more detailed theo retical treatment, we refer to lEmbrechts and Maejimal ( 20021 ) or 
Samorodnitsky and Taqqul (11994J ) . 

On the other hand, recently, extensions of the classical OU process have been suggested 
mainly on demand of applications. The generalized Ornstein-Uhlenbeck process given below is 
one with Levy processes and plays important role in Economics (pricing of Asian options , perpe - 
tuities and risk theory). For its theories and applicatio ns, we refer to e.g. ICarmona et all ( 2001 ). 
Erickson and Mailer! (120051 ). iLindner and Mailer! (120051) and Kliippelberg and Kostadinoval (2008). 
A multivariate extension is also considered in iKondo et al.l (120061 ) and lEndo and Matsuil (120081 ). 
Let {(£t, r] t )} t > be a bivariate Levy process and Vo be an independent initial random variable. A 
generalized Ornstein-Uhlenbeck (GOU) process is defined as 
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The stationarity and the convergence or divergence property have been intensively studied. If 
{£t}t>o an d {vt}t>o are independent and J Q e~^ 3 -dr] s converges a.s. as t — > oo to a finite random 
variable, V := {Vt}t>o has the stationary version (see e.g. Remark 2.2 of 
( 2005 )). The short memory property of V was also shown in Section 4 of 



Lindner and Mailer 



Lindner and Mailer 



(120051 ). 

Another extension of the original OU process is the fractional Ornstein-Uhlenbeck process, 
where FBM is used as integrator. An advantage of using the process is to realize stationary long 
range dependent processes. Let A > and an initial random variable Xq Gl 1 . A fractional 
Ornstein-Uhlenbeck (FOU) process is defined as 
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Here we need a non-semimartingale approach to cons t ruct s t ochastic integrals with FBM. We 



can find several u s eful t h eoretical tools in e.g. , iLinl (119951 ). iMikosch and Norvaisa! (120001 ) or 
Pipiras and Taqqul (120001 ) . ICheridito et al.l (120031 ) has shown that the FOU process is the unique 
continuous solution of a Langevin equation: X t = Xq — A j Q X s ds + Bf , t > and investigated 
its dependence properties. The main purpose of this paper is to construct a version of the GOU 
process which allows for long memory modeling by the use of a FBM. 



In order to define a generalized Ornstein-Uhlenbeck process we define a two-sided Levy process 

as 

t - i il if * " ° (3) 

where {Ct }«>o an d {£f }t>o are independent copies of {£,t}t>o- We work throughout with a bivariate 
complete probability space 

(ft := ft x x fi 2 , ^ := #x ® ^ 2 , P := Pi ® P 2 ) . (4) 

Let {£t}teR defined on (fti,«^i,Pi) be a Levy process and a FBM {P^} t6 K with Hurst index 
H G (0,1) defined on (^2,^2, P2) which is independent of {£t}ieR. A generalized fractional 
Ornstein-Uhlenbeck (GFOU) process with initial value Y G L 1 (ft) is defined as 

Y t : = e" 5t (y + J e?-dBf\ , t > 0. (5) 

If the initial variable satisfies 

Y = f e?-dB?, (6) 

J —00 

then, for convenience, we sometimes replace Y = {Y t } t > with 

Y t :=e-^ [ e^-dBf. (7) 



00 



The process Y is regarded as an extension of V given in ([T]) where the stochastic process of 
integration {r] t } t > is replaced with a {B]*} t > with PT G (0, 1) and also is regarded as an extended 
version of the FOU process where the integrand is replaced by the exponential of an independent 
Levy process £*. We should remark that Y has jumps caused by the process e~&. 

The paper is organized as follows. In Section [27T1 we recall the definition of Levy processes and 
summarize properties needed. In Sections 12.2.11 and 12.2.21 we review Riemann-Stieltjes integrals 
for functions with bounded p- variation and the stochastic integral in the L 2 (f2)-sense respectively. 
In Section [3] we investigate the existence of the integral in the GFOU process in order to justify 
the definition of the GFOU process. The stationarity condition and the second order behavior of 
the GFOU process are discussed in Section m and we observe the long memory property. Here 
we also examine stochastic integrals constructed by a single FBM, where £ in the process Y is 
replaced with B H used as the integrator. In Section [S] we obtain a stochastic differential equation, 
whose solution is given in form of the GFOU process. 

We use the following notations throughout. Write =' if equality holds almost surely. We will 
take the expectations for a bivariate process {{Z}, Z t 2 )} ie K. If the expectations only for a process 
{Zl} te $i is considered, we write its expectation as E z \. 

2 Preliminaries 

2.1 Levy processes 

In this subsection we introduce the setup for the Levy process. Let £ := {£t}i>o be a Levy process 
on M. with (a^, v^, 7^) generating triplet, where a^ > and 7^ G M are constants and a measure v^ 



(1 A \x\ 2 )u^(dx) < oo. 



onl\ {0} satisfies 



We call Ug the Levy measure of £. Then, the characteristic function of £ t at time t — 1 is written 
as 



exp 



at 2 

—z +i^z + 
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(e izx - 1 - izxl { \ x \< 1} ) vtr(dx) 
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For more on Levy processes and their properties, we refer to ISatol (119991 ). In later sections we 
consider several examples related the a-stable Levy motion with index < a < 2, denoted by 
£, a '■= {£f }t>o- K is a Levy process and its generating triplet is (0, u a , j a ) where 



u a (dx) :- 



C\X a dx on (0, oo) 



C 2 \x\ 



-l-a 



dx on 



-oo,0) 



with c\ > 0, C2 > and c\ + c 2 > 0. 

In order to define the in (jHJ), the variation of the Levy process plays an impor tant role. We 
give a brief summary based mainly on Section 5.4 of iDudley and Norvaisal (119981 ) and p. 408 of 
Mikosch and Norvaisal (2000). Define the p- variation for < p < oo of a process X := {X t } te ^ on 
[£i, £2] for ti < £2 in K as 



v p (X):=v p (X,[t 1 ,t 2 ]) 



sup 

A 



El^ 



X. 



(9) 



where A is a partition t\ — sq < S\ < ■ ■ ■ < s n — t 2 ol [ti, t 2 ] and n > 1. If v p (X, [ti, t 2 ]) < 00 for 
all t\ < t 2 , we say X has bounded p- variation, and if Vi(X, [ti,t 2 ]) < 00 we sa y it is of bound ed 
variation. Since every Levy process is a semimartingale v p (£) < 00 for p > 2 (see lLepingld (119761 )). 
We will state three useful results which characterize p-variation of Levy process in terms of the 
Levy measure. Unfortunately we can not find a result which uniformly characterize the variation 
in terms of the Levy measure. Assumptions and results are some what d ifferen t from paper to 
paper. The first one is a well-known result (e.g. Theorem 21,9 (i) of ISatol (119991 )). 



1. Bounded variation 

A Levy process {£t}i>o is of bounded variation if and only if 



(lc 



and 



[1 A \x\)v^{dx) 



< 00. 



The next res ult is a combina tion of Theorems 4.1 and 4.2 of iBlumenthal and Getoorl (ll96ll ) and 
Theorem 2 of iMonrod (119721 ). 



2. (a) p- variation of Levy processes 

Define j3 = inf{a > : J< i <x \x\ a v^{dx) < 00}. We call j3 the Blumenthal and Getoor 
index. If the Levy process {(,t}t>o has no Gaussian component (a^ = 0), then 



v p (£; [0, 1]) < 00 a.s. if p > j3 
v p (£; [0, 1]) = 00 a.s. if p < (3. 



(10) 



The result by iBretagnolld (119721 ) is a sharpened version of 2. (a) but with zero mean assumption. 

2.(b) p-variation of Levy processes 

Let 1 < p < 2 and {(,t}t>o be a Levy process without Gaussian component (a^ = 0). 
Then v p (£; [0, 1]) < oo a.s. if and only if 

(1 A \x\ p )^(dx) < oo. 



Otherwise v p (£; [0, 1]) = oo a.s. 



In particular for a-stable Levy p rocesses we have the result b y iFristedt and Taylorl (119731 ) which 
was stated in convenient form in lMikosch and Norvaisal (120001 ) . 



3. p-variation of ct-stable Levy processes 

Let {£"}t>o be a-stable Levy motion. Assume that 7 Q = for a < 1 and that the 
Levy measure is symmetric for a = 1. Then v p (C, a ) is finite or infinite with probability 
1 according as p > a or p < a. 

For the existence of the infinite interval integral in {Y t } t > given in (J7|) w e further need the 
behav ior of {(,t}t>o as t — > oo. Our assumption is that lim^oo^ =' +oo. iDoney and Mailer 
(120021 ) [Theorem 4.4] have obtained an equivalent condition of this in terms of the Levy measure 
z/£. Since se veral papers well exp l ain eq uivalent conditions, (see p. 72 of lErickson and Mallerl (120051 ) 
or p. 1704 of iLindner and Mallerl (120051 )) we do not mention it. Actually if lim^oo £ t =' +oo holds, 
we can assert a stronger result, which is more useful for our aim. 



Lemma 2.1 Suppose lim^oo^ 
t = to(uJi) < oo such that 



l-oo. Then for almost all io\ G Q± there exist 5 > and 
£t > St for t > t . 



The proof is a combination of Theorem 4.3 and 4.4 i nlDonev and Mallerl (20021). Concerning the 



integral J Q e ^ s ~dr] s in the GFOU process given by (Til) . lErickson and Mallerl (120051 ) have character- 
ized the convergence of improper integral j e~^ s ~dr] s in terms of the Levy measure of {(£*, Vt)}t>o, 
in which the condition lim^oo £ 4 =' +oo was used. 



2.2 Integrals with respect to functions with unbounded variation 

2.2.1 Riemann-Stieltjes integrals with p- variation 

We review several useful definitions of integrals of functions which have un bounded variation but 
bound ed p-variation. The excellent introduction to this area is given by iDudley and Norvaisa 
(119981 ). Let / and g be real functions on [a, b]. Define k = {u±, . . . ,u n } to be an intermediate 
partition of A = [a — s < Si < ■ ■ ■ < s n — b] given as in ([ 
i — 1, . . . , n. A Riemann-Stieltjes sum is defined as 



namely, s 3 -_i < Ui < Sj for 



Sns(f,g,A,K) 



1 1 

E 



f(ui) \g(si) - g(si-i)] 



Then we say that the Riemann Stieltjes integral exists and equals to /, if for every e > 0, there 
exists 5 > such that 

\Siis(f,g,A,K) -I\<e 

for all partition A with mesh ma x (s, : — s ,;_i) < 5 and for all intermediate partitions k of A. The 



r 

follow ing; theorem was pro yed bvlYouna (119361). (See a lso Theorem 2.4 of iMikosch and Norvaisa 
( l2000h or Theorem 4.26 of iDudlev and Norvaisal dl998h .) 



Theorem 2.1 Assume f has bounded p '-variation and g has bounded q-variation on [a, b] for some 
p, q > with p~ x + q~ l > 1. Then the integral J fdg exists in the Riemann- Stieltjes sense, and 
the inequality 



fdg 



holds with C 
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< c p , q (Mf)) 1/p (v q (g)) 1/q 
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2.2.2 Integral with respect to FBM in L 2 -sense 

Another definition of the integral is in L 2 (f2 2 )-sense. Stochastic integrals wit h re s pect t o FBM is 
sometimes defined as the L 2 (f2 2 )-limits of integrals of step functions (see e.g. iLinl ( 199a )). We see 
this when B H is the Brownian motion B 1 / 2 . If a function f(u) -."SL^ 1 ^ satisfies f(u) G L 2 (R), 
there exist step functions 

n 

fn{u) := ^ fi l {u,<u<u t+1 ], -OO < Ui < . . . < U n+1 < 00, fi 6 R, U E N. 

t=l 

such that f n converges to / in L 2 (R). Then the integral J R f(u)dB u is the L 2 (f2 2 )-limit of the 
integrals of step functions 

» n 

Uu)dBl' 2 = Y,n{B l Jl-BU 2 

Jr i=1 v 

since L 2 (M) and L 2 (Q) are isometric and their inner products are equal, namely, 



E 



(f n (u) - f{u))dB]/ 2 



(fn(u) - f{u)fdu. 



Pipiras and Taqqul (J2000l ) have investigated a similar characterization for the integral of B H 
when H ^ |. We apply this to the existence of the improper integral in the GFOU process {Y t } 
afterward. Define the linear space 



Sp{B 



H\ 



A : / f n (u)dB u — ► X in L (fi) for some sequence (f n )neN (step functions) 



Pipiras and Taqqul (120001 ) have analyzed the functional space of the integrand f(u) in which it can 
be asymptotically approximated by f n (u) and J R f(u)dB^ is well-defined. For H G (0, s) they 
succeeded in specified a Hilbert space of functions on the real line which is isometric to Sp(B H ). 
However, for H G (~, 1) they had difficulty in finding the corresponding isometric space, and 



as second best they analyzed inner product spaces in which the integral with B H (H G (1,1)) 
is well-defined. We give only one such inner product space and its inner product for B H with 
H G (s, 1)- Other inner product spaces do not seem to work for our purpose since they require 
e.g. characteristic function o f f or fractional deriyative of / which do not exists in our case where 
f(u) = e^ u . (See Section 7 of iPipiras and Taqqul (J2000l ).) The space is 



|A|*={/:<|/|,|/|)| AF <oo}, Hef^lY (11) 

where 

(f,g)\A\ H = c H / / f(u)g(v)\u-v\ 2H ~ 2 dudv 

Jr Jr 

is the inner product with ch '■= H(2H — 1). 

3 Existence of the integral 

In this Section we analyze the existence of the integral in the GFOU process given in (j5j). The 
definition of the GFOU process includes an integral with respect to FBM. Since paths of FBM are 
of infinite variation and since FBM with H ^ | is not a semimartingale, the stochastic integral 
with respect to FBM (H ^ |) is not an Ito integral. Additionally, the integrand of the GFOU 
process is random and the infinite interval integral ([6]) is needed for its stationarity. We apply two 
approaches of the integral in Section 12.21 in order to cope with these problems. 

Proposition 3.1 Let B := {B^} te ^ be a FBM with H G (0, 1) and £ := {£t}ieR be an indepen- 
dent two-sided Levy process. Assume that {C,t}teR has bounded p-variation for < p < oo. Then 
the integral J e^-dBf , < t < oo exists in the path-wise Riemann-Stieltjes sense if - + H > 1. 
Furthermore with q > jj and C P:Q = ^ n>1 n'^ 3 +q > we have 



t 

e^-dBf 
o 



< C P J sup <A (v p (£)) 1/P (v q (B H )) 1/q P - a.s. 
\se[o,t] J 



. a€[0,t] 

The proof of Proposition 13.11 is obvious from the continuity of the exponential function and we 
omit it. 

Remark 3.1 If p = 1 in Proposition 13.11 £ has bounded variation and we can define path- wise 
integrals for all B H with H G (0, 1). On the other hand if H > |, we can define the path- wise 
integral for any integrands £ since v p (£) < oo for p > 2. 

The reason why we need the path- wise definition besides the L 2 (fi)p-approach is that we want 
to define the integral for H G (0, |), which we could not obtain in the L 2 (f2)-approach. It also 
gives several useful tools easily, such as integration by parts or chain rule for analyzing stochastic 
differential equations in Section El 

Example 3.1 Set £ := £ Q be an a-stable Levy motion with index a defined in Subsection 12.11 
Then under assumption in 3. p-variation of a-stable Levy processes and the assumption 
- + H > 1. the integral e~^ f e*"-dB^ exists in the path- wise sense. Hence the GFOU process 
is well-defined. Both £ a and B H have infinite variation and the former is known to be extremely 
heavy tailed process. 



Recall that we write Y t = Y t (u)i,oj2) to emphasis that it is a function from probability space 
fllj) and that we write Y t if initial value Yq satisfies ([6]). Define approximating step functions of 
the integrand e~^ t+ ^ u " of Y t as 

n 

ft,n(u]^l) '■= 2 J fi( u 'i u; l)^{uf_,<u<uV'}, 
i=l 

where 

Here w"'s are points in [— N n ,t] such that — N n — Uq < u™ < ■ ■■ < u™ — t and as n — > oo 
max (u™ — u™^) j and iV n " oo. By integrating / 4j „ with respect to Bf, we also define 
approximating sequence as 

Z?{u l7 u 2 ) := I / tl „(u;wi)d£*(u; 2 ) 

>/ — oo 

= J]/ i ( W ;a; 1 )(^ 2) - J B5_ i (a; 2 )). 

In the following theorem we define the integral J_ e~^~^ u - dB^ as the limit in probability of the 
Z 4 n as n —* oo. The reason why we need this approach is that with only path- wise definitions we 
find difficulty to treat improper integrals. For the existence of the improper integral we should 
consider long time (t —> oo) behavior of both {B^} teK and {£t}t g R path-wisely, which seems to 
be not an easy task. Additionally, this approach is well-matched with the analysis of the second 
order behavior. 

Theorem 3.1 Let {£>^} teK be a FBM with H 6 (|, 1) and {£,t}tm be an independent two-sided 
Levy process. Iflimt-^oo^t = +oo, then for each t > Z™ given in (dp converges in probability 
to a limit defined by Y t and which does not depend on the sequence u™. If further E[e~ 2 ^] < 1, 
then Z™ converges to Y t in L 2 (Q) and it follows that 



e -&-Z«-)dB* 







E[Y t ] = E 
and for < s < t, 

E[Y s Y t ] = I I E^[e-^ + ^--^+^-]\u-v\ 2H - 2 dudv. (12) 

J —oo J —oo 

Remark 3.2 (a) For H e (0, |) we could not validate the existence of the improper integral in 

Y t . 

(b) Under the only condition that lim^oo^ =' +oo, we cannot prove the L 2 (Q) convergence. 



c) The condition E 



Lindner and Mallerl (120051 ) 



< 1 implies lim^oo^t =' +oo, which is shown in Proposition 4.1 of 



Proof of Theorem \3.1\ 

We check that for each u% the integrand e~^ Ul ' + ^ u -^ 1 ' E \A\ H given in (iTTj) . Since ^t(^i) is 

constant, we drop it and only show e^" - ^ 1 ^ G |A| H . The function e^ - ^ is caglad and bounded 



on any finite interval. Due to Lemma 12.11 and to the symmetry of two-sided Levy processes there 
exists T{uji) < such that for all u < T(ui), £ u _ < 5u where 5 is some positive constant. Then 
for each ui, 



t r-t 



e €„- («i)+€„- («i) I u _ v \^-2 dudy < ^ 



OO J — OO 



is obvious. Hence we can utilize L 2 (f22)-integral theory in Section [2.2.21 Namely, for each t > 
and for each fixed u>i, Z™(ui, •) converges in L 2 (f22, i^)- Moreover Z™ converges in probability on 
(fl, P) for each t > since sequence Z" satisfies the Cauchy criterion, as seen by 



n,m—* 00 



lim P(|Z t n -Z t m | >e) 

00 

lim / / l { \ Z r, _ zrl>e} dP 2 dP 1 
[ lim P 2 (|^> 1 )-ZrM|>e)cLP 1 = 0. 



The limit is called Y t and it is J^i ® ^ 2 measurable for each t. 

Now, with E[e~ 2 ^] < 1, we prove the L 2 (f2)-convergence. We have Pfe - ^ 1 ] < 1 as well, hence 

9 1 := -\ogE[e-^] >0 and 6 2 := -logP[e" 2€l ] > 0. 

Then using the covariance of the FBM, we have 



P[(Z 



D 2 ] 



E 



-2£t+£ u n , +£„ 



r rf f rf I z — 1 j j — 1 



n n 



t=l 3=1 



e '- 1 -J- 1 



( - |< - u™| 2// + \u n 



+\<-^-i\ 2H -K-i- u U\ 2H )} 



i-1 u j I 



ch 



t ft 



T n ( Ul v)\u-v\ 2H ~ 2 dvdw 



(13) 



00 j —00 



where 



!?(«,«) 



n n ^ 

Z^z_^ 1 1 K_i<«<«" 1 «"_i<^<«" 1 i>i} e " 
i=i 3=1 ^ 



■«a(t-«?_i)-ffi(«?_i-«7_i) 



"t" -*-{«"_! <«<«",«?_! <v<u" i<j} ( 



-fe(t-«J.iK(«J.i-«JLi) 



which obviously converges point-wise to 

1 r ^ n(9 -e 2 (t-«)-ei(«-w) , 1 -0 2 (t-i>)-<Mv-u) 



Since |r"(v, w)\ < M' for some M' > uniformly in n, t < 0, (TT3l) is bounded by cnM't H . 
Furthermore according to usual Fubini's theorem, it also follows that 



E [(Y t f] 



c H E^ 



(Y t (u 1 ,uj 2 )) 2 dP 2 dP 1 



e -2ttHu-Hv- \ u _ v \ 2H - 2 dudv 



■2c 



u 



oo J — oo 
t rt 



-e 2{ t-u)-e l{ v-u) l{u ^ u _ v \2H-2 dudv < ^ 



oo J — oo 



Observe that this integral i s finite. Acco rdingly E[(Z 



n\2l 



E [(y t ) 2 ] as ri -> oo. Now we 
can apply Theorem 4.5.4 of IChung! (1200 ll ) to Z™ and obtain the L 2 (r2)-convergence. In conse- 
quence i?[y s y t ] turns out to be finite and equation (fT2"j) follows from Fubini's theorem. Finally 
E[Z™] =0, fiGN implies E\Y t] = 0. Hence the proof is complete. □ 



The process {Y t }teR obtained in Theorem 13. II is the GFOU process with initial value Yq. We 
close this section with the following concluding Remarks. 

Remark 3.3 (a) In both Proposition 13. II and Theorem 13. II £ is independent of B H and we have 



e~^ / e^-dBf? 



e b+€»-dB**, -oo < a < b < oo. 



< oo a.s. Q. Therefore 



This is not allowed in usual theory of stochastic integrals related to semimartingale (IProtter 
( 120041 )) since {£*}teR is not adapted. 

(b) From Theorem 13.11 linu^,^ £+ =' +oo implies |lo| = /_ e^ B ~dB^ 

together with Proposition 13. II we can also treat \Y t }t> n = [Y t + Y } t > path-wisely 

(c) We should mention that lErickson and Mallerl (120051 ) [p. 81] gave another idea for improper 
integrals of a function of Levy processes with respect to FBM: J °° g(£ t )dB^ . Investigation of their 
idea is also our next concern. 



4 Stationarity and Second order behavior of GFOU pro- 
cesses 

Here we investigate the strict stationarity and the second order behavior of the GFOU process 
Y := {Y t } t > . Since we could not validate the existence of Y with Hurst parameter H £ (0, |) 
and since our main concern in this paper is the long memory case, we confine our results to the 
case H £ (|, 1) throughout this section. 
The stationarity of Y is as follows. 



Proposition 4.1 If Hindoo £t 
strictly stationary. 



+oo, then Y t exists for all t and the process Y := {Y t } t >o is 



Proof of Proposition \4-l 



Let < tx < t 2 < . . . < t m < oo, m £ N and h > 0. We use the sequence Z™ given in ( 1T21) . Since 
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both {£j}t e R and {B^} t& R have stationary increments, does the pair {(£ t , £? t H )} ie R as well because 
of independence. Thus, 



3=1 
n 

• -tt+h+iv" ,+h ./ tdH tdH \ 

e j-i (B u n +h - B u n +h ) 



= 2*_ 1+h , l<i<m, 
simultaneously in z and we have 

l-^ii) • • • > ^t m ) — y^ti+hi ■■■■> Zj t m +h)- 

Then by virtue of Theorem 13.11 as n —* 00, (ZP , . . . , Z t n ) converges in probability to (F tl , . . . , Y" in ) 
and (Z? 1+h , ..., Z? m+h ) converges in probability to (Y tl+h , ..., Y tn+h ). This yields 

\Y tl , • • • 5 ^t„) = (Yti+h, • • ■ ) ^t„+/i) 
and the conclusion holds. □ 

Remark 4.1 (a) The logic in the proof works even in the case H e (0, |). If the assumptions of 
Proposition 13. II are satisfied and if the integral Yq = J^e^^dB^ 1 with H e (0, |) exists, then F 
with He (0, 5) is defined and is strictly stationary. 

(b) In connecti on with the GOU pro cess {V}t>o given in (JTJ). Proposition 14.11 corresponds to 
Theorem 3.1 of ICarmona et al.l (l200lh . where £ and 77 are independent, and under conditions of 



a.s. convergence of the integral L e ^ a ~dr] s and lim^oo^t =' +00 the stationary version exists 
and equals in distribution to \4o- 

Next we investigate the second order behavior of Y and derive the auto-covariance function 
explicitly. What should be remarked is that while the aut o-covariance function of the GOU process 
V given in ([T]) decreases exponentially (Theorem 4.2 of iLindner and Mallerl ( 20051 )). that of the 



FOU process \X^\ t >n given in (TTJ) decays like a power function (Theorem 2.4 and Corollary 



2.5 in ICheridito et al.l (120031 )). Since Y is regarded as a version of GOU processes and FOU 
processes, this investigation is interesting. We utilize results in Theorem 3.1 and obtain Theorem 
14. II and Corollary 14. 1 1 below. Note that even the existence of Y and the equation (Tl2|) are obtained 
several difficulties still lay in calculating the auto-covariance function. The integrand in (1121) is 
regarded as exponential moment of of 4 dependent random variable, i.e. E^[e~^ s ~^ u - + ^~^ v -'] and 
dependencies of these variables are different in the order of s, u, t and v. We also require that 
after E^[e~^ 3 ~^ u ~ + ^~^ v ~^] is calculated the double integral in ( TT2l) has a suitable representation 
for our purpose. In the proofs of Theorem 14.11 and Corollary 14.11 we will see how to get over these 
difficulties. Proofs of Theorem 14.11 and Corollary 14.11 are given in Appendix |X] since they require 
a lot of technical and tedious calculations. 

Theorem 4.1 Let {Bf} te ^ be a FBM with H e (|, 1) and {£,t}tm. be an independent two-sided 
Levy process. Suppose that E[e~ 2 ^] < 1. Then the stationary version Y := {Y t }t>o exists and for 
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any s > 0, t > 0, we have 



Cov(Y t ,Y t+s ) = 
2e~ 6lS 



Cov(Y Q ,Y s 



CH 



9 2 9 1 



2H- 



r r(2H-i; 



29{ 



211 



T(2H- 1) 



is s 2H-l 



29 lS 



29 1 (2H- 1 

2e~ 0lS 



iFtfH -1,2H; 9is) + t^2h T ( 2H ~ *> Ois 



201 



ch 



ns 



9 2 9 x 



2H- 



T T(2H 



291 



211 



Y(2H 



-,2H-l 



OO 71+1 



iV ' n=0 fc=l 



?lS )-(n + l)_ ( _ M -(n +1 )7(n+l,M 



where 9± := — log(-E[e - ^]) > an d 6?. := — logf-E^e" 2 ^ 1 ]) > 0. i/ ere 7(-, •) and T(-,-) are incom- 
plete gamma functions in 8.350 of lGradshteun and RuzhiA i 200(a) and iFi(-, •; •) zs t/ie confluent 
hyper- geometric function in 9.210 of lGradshteun and Ruzhitt \2000\ ). 

Note that since 7(77, + 1, s#i) — > r(n + 1) = n\ as s — > 00, we have 

00 2n— 1 



n=l fc=l 
2H-2\ 



= 0(s 

This conclude that F with H G (|, 1) is a long memory process. While we obtained Cov(Y t , Y t + S ) 
using special functions in Theorem 14. 1[ it is mainly for numerical purpose since for such functions 
useful softwares are available. 

Next we investigate long time dependence of {Y t } t > with the initial value Y := X G L 2 (Q) 
where X is independent of {£t}t>o and {B^} t >Q. 

Corollary 4.1 Let Y := {Y t } t > , H G (|, 1) be a GFOU process with the initial value X G L 2 {VL), 
where X is independent of £ := {Ct}t>o and B H := {B^} t > . Then for fixed t > as s —» 00. 



00 2n— 1 



Cov(Y t ,Y t+s ) = Hj2]\(2H-k)9^ n {s 



2H-2n 



hi 



(t 



\2H-2n 



} + 0(s 



2H~2N-2\ 



0(s 



n=l k=l 
2H-2\ 



We see what happens to the second order behavior of the process Y if £ in Y is replaced with 
B H which is the same process as the variable of integration. Although we expect long memory 
this does not hold. Note that we need only the probability space (f2 2 , J^2 5 -F2) here. For H G (|, 1) 
and the initial random variable X G L 2 (Q 2 ) independent of {B^} t( zj&, define 



W t 



-B? 



X+ / e B "-dB* 



To analyze {Wt}t>o we use the path-wise integral theory (see Subsection 12.2.11) . Let / be contin- 
uous different iable and F{x) = F(0) + J* f(y)dy. Then with H G (|, 1) it follows that 



F{B? ) - F{B* ) = f f{B*)dB« 

Jo 



a.s. 
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By setting / = e 4 in above we obtain 

W t = l + e- B * B (X-l) 



Proposition 4.2 Let H E (~, 1) and t, s > 0. Define M x := (E[X) - l) 2 and M 2 : = E[(X - V 
The process {Wt}t>o has the following auto-covariance and correlation functions. 

Cov(w t ,w t+s ) = e H* 2i? +(*+^}| M2e |{^+(t+ S )--^}_ Ml } 

M 2 e^ 2H ~^-M ie -^) 2H 



Covv(W t ,W t+s ) 



VM 2 e tH - M 1 y/M 2 - M 1 e-(*+ s ) 2 
0(e _ 2 s ) as s — ► oo. 



We also consider the drift added process 



W t = e^ B " + ^ (x + f' e B " +au d(B B 



where a > and X G L 2 (Q 2 ) is independent of {Bf}t>o- Even if a drift is added, usual path-wise 
integral works and 

W t = l + (X- i) e -(B?+*t) 

holds. Then the auto-covariance and correlation functions of {H^t}t>o are calculated is a similar 
manner and become 

Cov(W u W t+s ) = e- a ^Cov(W t ,W t+s ), 
Cott(W u W t+S ) = Coiv(W u W t+s ). 

Thus our conclusion here is that even if a drift is added it does not have long memory. 

5 Stochastic differential equation related with GFOU pro- 
cesses 

We analyze a stochastic differential equation of which a solution is given by the GFOU process. Let 
U := {U t } t > be a Levy process with generating triplet (ay, vy, 7c/)- Assume that the Levy measure 
vy has no mass on (— oo, — 1]. The Doleans-Dade exponential of U t is written as £{U t ) = e~& 
where 

6 = -u t + ^ t - J2 (Mi + *u.) - AU S ) . 

0<s<t 



See Section 2.2 of lErickson and Mallerl (120051 ). Here £ f is the Levy processes. 



13 



Proposition 5.1 Under the assumption in Proposition ^ '. 11 GFOU {Y t }t>o with the initial value 
Yq G L x (Q) satisfies the stochastic differential equation; 



dY t = Y t ^dU t + dB? 



t ■ 



(14) 



where £{U t 



Proof of Proposition I5.il 

Since the condition of Theorem 12 .11 is satisfied, the integral f Q B^_de^ s also exists in the Riemann- 

Stieltjes path-wise sense. We use the integration by parts formula to Y t and obtain 

Y t = e~^ [Yn - [ B H deH + B?. 



If we put Qt '■— Y t — Bf, the equation above becomes 

Qt = e~ Ct (qo - J B?_d£ 

where Qo = Y . Since {e^} t >o is a semimartingale and {B^} t >o is continuous and adapted, the 
process {Qt}t>o is also semimartingale. We set R t := e~^ and S t : = Q n — J Q B f de^ a . Then the 
integration by parts formula for semimartingales (e.g. Corollary 2, II of IProtterl (120041 )) yields 

Qt ^ Qo — RtSt — RqSq 

ft rt 

R s -dS s + / S s ~dR s + [R, S]t — RqSq 



- f e"*-.B 5 _de & + J Q s ^ s -de^ s - f B^d[ e -^e^} 
Jo+ Jo+ Jo+ 

[ Q s ^dU s - [ £JL(e-*-de*'-d[e-*,e*] 8 ). 

Jo+ Jo+ 



Observe the relation between e & and U t ; 

1 



e 6 e 6 



/ e is ~de ?a + / e ia ~de is + [e*, e € ] 
Jo+ Jo+ 



dU„ 



n+ 



-de is 



[e^e-%. 



Using this we obtain 



which is equivalent to 



The proof is now complete. 



Qt - Qo = I (Q s - + B?_) dU s 

J0+ 



Y t -Y 



Y_dU s + B 



it 



D 
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Remark 5.1 The Levy measure of {£,t}t>o is obtained from that of {Ut}t>o] 

u^((x,oo)) = u u ((—oo,e~ x — 1)) and z^((— oo, — x)) = vu((e x — 1, oo)). (15) 

See again Section 2.2 of lErickson and Mallerl ( 20051 ). Hence if vu is concretely given, using criterion 



of p- variation in Section 12.11 we can check the condition of Proposition 13.11 

The following technical Lemma is not difficult but useful for the existence of {Y t } t > Q which is 
directly constructed from the stochastic differential equation (JHJ). The poof is only a calculation 
and we omit it. 

Lemma 5.1 Assume that {Ut}t>o is a Levy process and that {£,t}t>o satisfies £(Ut) = e~& . Then 
for < 5 < 2, convergence and divergence of 



\x\ u^(dx) and / \x\ uu(dx) 

x\<l J\x\<l 



are equivalent. 



Example 5.1 As an example we consider the stochastic differential equation ( Tl4|) . where Ut is 
given by an a-stable Levy motion £° (see Section |2~T|) . From remark above the Levy measure z/g 
is given by 



v ( dx \ = f c 2 (l-e x ) ^ a e x dx on (0, oo) 

\ c\ (e~ x — 1) a e~ x dx on (— oo,0). 
Observe that 

u^(dx) ~ \x\~ 1 ~ a dx as |x| { 

and hence variation property of £ t is the same as that of Ut = £"■ As a result f p (£) is finite if 
p > a. 



A Proofs of Section 3 



Proof of Proposition^. 1 



From the stationary version Y is definable. By virtue of the stationarity of Y and Fubini's theorem, 
we have 

Cov(Y t ,Y t+s ) = Cov(F f ,F t+s ) 

"0 pO 



c H f f Es[e-&-*"— ^- ] }\u-v\ 2H - 2 dudv 

J — oo J —oo 

+c H [ [' Et[e-&-^—^]\u-v\ 2H - 2 dudv 

J -oo JO 



: = I + 11. 

First we consider the integral I. The independent increments property of {£t}t e K gives 
£;r e -(£s-£u--&>-)i — E\ e -{£s-£u-£v)i 

= E[l {v > v} e-&-M- 2 to-*J-&-*J] 

+E[l {u<v} e-^o)~2(t;o^)-(^u)} 

= l {u > v} E[e-^-^\E[e- 2 ^-^-^\E[e-^^-^\ 

+ l {u<v} E[e^^ o) }E[e~ 2 ^-^^ o) }E[e~^~^°\ 

15 



The integrand E[e ^ s ^ u ~ €■"->] is symmetric with respect to u and v, and hence 

-0 r 



\ = 2c H e^ s f I l {u > v} e^ u+e ^-^\u-v\ 2H - 2 dudv. 

J — oo J — oo 



Then further calculation shows 

»0 /-O 



I = 2c H e^ s f I l {u > v} e 9 * u+9 ^-^\u-v\ 2H - 2 dudv 

J — oo J —oo 



-oo ./ — oo 

(By change of variables; x — u — v) 

-0 poo 



/U ^oo 

e d2U du / e" eiX x 2H " 2 rfx (16) 

■oo J0 

= ZCHe T(2H-1). 

Next we consider II. A similar conclusion as above gives 

j5;r e -(k-6'-)-(£o-&.-)l _ £r e -(&-&0-(&-6u)| 

= e - 6 'i( s - 1 ')+ 6 'i («+ u ) 

_ e -0is+0i(it+t;) 

and we have 

II = Cff e" eiS / e eiM /" (u - v) 2H ~ 2 dudv (17) 

./-co Jo 

(By change of variables; x = v — u) 
= c H e~ 9lS I e 9lU I U e 9liu+x) x 2H ~ 2 dudx (18) 

J —oo J s—u 

poo r0 

= c H e~^ / l { _ x<u<s _ x] e 2e ^*x 2H - 2 dx 

Jo J -oo 

(PS /*0 /*oo /* s — X 

/ e e ^x 2H ~ 2 dx I e 2e ^du + / e 6lx x 2H - 2 dx / e 2 * 1 ^ 
Jo J— a; Js Jx 

= CHe -^ (-L f e ^ x x 2H - 2 dx - -L /"°° e" 91 V H " 2 rfx (19) 

V 2 ^i io 2 ^i io 

D 2#is /-oo 

4 



20is /-oo \ 



s 2H+1 

2e 1 W^T) lFl{2H - h2H;6ls) -W 



c He~ 0lS ( 0/Worj ^ iFi(2ff - 1, 2#; l5 ) - -^r(2# - 1] 



e 2e lS 



T(2H-l,9 lS ] 



29 2H 
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Combining I and II, we obtain the first assertion. 

The series representation of the incomplete Gamma function, i.e. 



T(a,x) = x 



( °° 1 

\ n=l 



[a — n) 



gives 



_%H-\ °° n + l 



T(2H - 1,9,.) = — - _ Y. II( 2g - *)(«")-'• 

n / n =0 fc=l 

We apply the binomial expansion 

(2H - 2){2H - 3) • • • (2H - n - 1 



+i) 



(20) 



2H-2 



(1-u/s) 



to the representation 

e -6is s 2H-l 

20i(i/-l) 



E 

n=0 



F 1 (2ff-l,2ff;0 1 3) 



in 



s. 



(0 < u < s) 



hs r s 



e e lXx 2H-2 dx 



(By change of variables; u = s — x) 



J2H-2 ps 



e~ yiu (l - u/sY H - 2 du. 



26*i J 
Then we exchange the infinite sum and integral by usual Fubuni's theorem and obtain 



e -e lSs 2H-i 
2GAH-1) 



iFi(2#-l,2#;#is) 



~.2H-1 °° n+l 



■—4 rV ~\(2H-k)(-sO 

2d x (2H-1 )^i-i- K A 

J-V / n =0 fc=l 



N -( W +1) 7(W+1,M 

lj j ■ 

n! 



(21) 



Thus substituting these expansions ( !20l) and (121J1 in the previous representation of covariance we 
obtain the result. □ 



Proof of Corollary \4-l\ 

Since X is independent it follows that 

E[Y t ]E[Y t+s ] = (E[X}) 2 e^ 2t+s \ 

We divide E[Y t ,Y t+s ] into piece as follows. 

E[Y t Y t+s ] 



(22) 



E 



Xe^ + j e-^-^dB^llxe-^ + f e^+'-^dBf 



E[X 2 e^ tHt+s) ] + E 
+E 



Xe 



-& 



t+s 



-(£t+a-tv-) flgH 



o 
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£ /-t+s 



LJO </0 



e -&-i«-+Z*--Zv-)dB*dB* 



The first term is calculated as 



E[X 2 e- { ^+M} = E[X 2 }E[e~^ t+ ^ t+2 ^] 
= E[X 2 ]E[e-^}E[e- 2 ^} 
= E{X 2 ]e~ 8lS - e2t . 



(23) 



From Theorem 14.11 the second and the fourth terms are 0. We only need the last term, which is 
calculated as 



Cov{Y t ,Y t+s ) - E 
-o 



ii 



t+s 



-E 



-E 



e -(6-e-) dB f / e -(^-e--)dBf 

oo J t 

t+s 

I ^~tu-) dB H I e -Ht+s-^-) dB H 



u 



V 



■oo J — oo 

e -^~tu-) dB H / e -^+s-i^) dB H 



Cov(Y t ,Y t+s ) - {e- 9lt Cov{Y ,Y t+s ) - e~ e ^ Gov (Y ,Y t )} 

-e- 9lS E [e-^YoYt] - e~ eiS E[e~^Y \Y t 



°1 

t h 



where Y,° is Y t with initial value 0. 



By adding up (1221), (T23T) and ()24jh we have 



Cov(Y u Y t+s ) 
= Cov(Y t ,Y t+s ) - e- 9li Cov(F ,F i+s ) 

+e~ eis ( E[X 2 ]e- 6lt - (E[X]) 2 e- 26lt + e - 6lt Cov (Y ,Y t 



-E[e-^Y Y t }-E[e-^Y Y t ° 



Hence via Theorem 14.11 the conclusion follows. 



(24) 



□ 
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